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Summary 


The research is concerned with developing a new approach to enhancing fault tolerance 
of flight control systems. The original motivation for fault-tolerant control comes from 
the need for safe operation of control elements (e.g. actuators) in the event of hardware 
failures in high reliability systems. One such example is modem space vehicle subjected 
to actuator/sensor impairments. A major task in flight control is to revise the control 
policy to balance impairment delectability and to achieve sufficient robustness. This 
involves careful selection of types and parameters of the controllers and the impairment 
detecting filters used. It also involves a decision, upon the identification of some failures, 
on whether and how a control reconfiguration should take place in order to maintain a 
certain system performance level. 

In this project new flight dynamic model under uncertain flight conditions is considered, 
in which the effects of both ramp and jump faults are reflected. Stabilization algorithms 
based on neural network and adaptive method are derived. The control algorithms are 
shown to be effective in dealing with uncertain dynamics due to external disturbances 
and unpredictable faults. The overall strategy is easy to set up and the computation 
involved is much less as compared with other strategies. Computer simulation software is 
developed. A serious of simulation studies have been conducted with varying flight 
conditions. 
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1. Introduction 

A growing importance has been given in the last few years to the design of a flight 
control system able to accomplish the reconfiguration of the aircraft following battle 
damage [l]-[5]. While there exists many control methods for flight control, most of the 
existing methods are based on "healthy” aircraft model in that the control algorithms are 
derived from system dynamics without any faults in the system. Evidently, automatic 
control plays a crucial role in achieving these features. 

This work is concerned with the development and analysis of neural network based 
control algorithms for flight vehicles with faulty dynamics. A neural-adaptive method to 
enhance flight vehicle stability under unpredictable faults is presented. The method is 
based on the dynamic model coupled with the effect of unpredictable faults and 
disturbances. It is shown that with the proposed strategy the effect of faults on the system 
performance can be effectively suppressed. The novelty of the proposed approach also 
lies in the fact that it is fairly easy to set up and the computation involved is much less as 
compared with other strategies. An example is used to verify the validity of the proposed 
approach. 


2. Problem Formulation 

The major thrust of this research is directed toward the stabilization of flight vehicle 
subject to subsystem failures. The dynamic model considered is of the following form, 

M(.)r+N(r,r)+d(t)=Bu (1) 


with 


M (.) = M o ± f$(t - Tj )M f (.) 


N(.) = N 0 ±/3(t-T f )N f (.) 


where M e is the inertia matrix depending on the mass property of the vehicle 
structure, re R n denotes the coordinates/states of the system, Ne R'”™ represents the 
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term related to damping and stiffness properties of the system, de R n represents the 

external disturbance, we R m denotes the control signals of the actuator, and Be R nxm is 
the control input influence matrix reflecting the location of the actuator distributed along 
the vehicle structure. In most cases B is a matrix with 0 and 1 as its elements. 

The notation “±” is used in the matrices M and N to emphasize the fact that both the 
direction and the magnitude of the uncertainties due to modeling approximation and 
possible faults are not available for control design. 

Note that the f(.) t and fi(t-Tf)Bf{.) are used to model the 

unknown portion of the dynamics due to a fault occurring at time t>Tf . The function 
takes the form 


0 t<T f 




which covers both jump and continuously time-varying (ramp) faults. For the system to 
admit a feasible solution under faulty condition, we need to assume that there exists a 
positive constant c 0 > o such that 

J Pit -T / )M / (.)|+c 0 <|M 0 | (2a) 

This assumption implies that the fault condition under consideration is within attackable 
range in that it does not dominate the value of M under normal conditions. In such a case, 
the symmetric and positive definite property of the mass matrix still holds, i.e., 

Af (.) = M (.) r > 0 (2b) 

Whenever a fault occurs, system dynamics may be changes in many different ways. In 
such a case, our primary objective is to design a proper control to maintain the stability of 
the system. Note that essentially two system modes are involved under any fault: critical 
modes (denoted by r c ) and non-critical modes (denoted by r ^. ). The critical modes 
should be stabilized explicitly for they are directly related to system safety. 


P(t-TA 
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For control design we assume that the rank of B is m, which implies physically that there 
are m actuators mounted at appropriate positions in the vehicle for fault accommodation. 
Obviously, only the mode, with which there is an actuator, can be actively damped out. 
Therefore, to make full use of the available actuators, we consider the case of dim( r c ) = 

m, thus dim( ) = n-m. Under the collocation condition (i.e., sensors and actuators are 
located at the same places), we have 

r c = B T re Q c c R m , r m = P jr e c R n ~ m (3) 

where Pj is a projection matrix. Clearly, £i. = Q c <BQ nc , here £2 denotes the whole 

coordinate space of practical interest. In this study, our objective is to develop a control 
scheme that reduces both |r c | and |r c | to an acceptable level, i.e., 

|r,»S £l , |K || S 

where e x and e 2 are small numbers. A strategy achieving this objective is developed in 
the next section. The following lemma is needed to develop the control strategy. 

Lemma: 

Let b > 0 be a constant chosen by the designer, and define 

s = (B T M~ l B)~ l B T (f+br ) (4a) 

If ran(B) = m , then 

(i) | ^| — > Oas t—>oo implies that |r c J -» Oand ||r c | — >0 as t — > <» . 

(ii) |H| <c<°° as t>T c implies that ||r c || < c r < «> and ||r c || <Cj.<°° as t— > «> . 

Proof: 

Let 

Z = r c +br c (4b) 

then from (3) and (4a)-(4b) we have 
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s = (B T M~ l B)~ l z 


(4c) 


With rank(fi)=m and the property of M as in (2b), it follows that b t m 'b is symmetric and 
positive definite. Therefore |s||— »0as t— implies that |z|| — > 0 as t— By the 
definition of z as in (4b), the results as stated in the lemma are readily established. 

As motivated by this lemma, we first construct the following signal 

z = r c +br c = B T (r+br) (5) 

where b> 0 is a positive number chosen arbitrarily, then we focus on making z go to 
zero or be less than a small constant which is determined by the required precision. 

Note that this signal is available as long as r c and r c are measurable. From (1) and (5) it 
can be shown that 

i = P(.)u + L(.) (6a) 

where 

P(.) = B T M~ l B (6b) 

U.) = bB T r-B T M'\N+d ) (6c) 

It should be emphasized that P and L are unavailable due to the fact that the structure 
contains fault and external disturbance. This calls for a more dedicated control design 
approach for such a system. 

It is interesting to note that the system model (1), although uncertain, possesses the 
following properties. 

Properties: 

1) M~ l =M~ T > 0. 

2) If B is full row rank, then there exists constants c min and Cmax such that 

c min M 2 ^ Z T B T M~ l Bz < C max | zf 
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T _i 

where c,^ > c min denotes the maximum and minimum eigenvalues of B M B 
respectively. 


Proof: 

The proposition can be easily shown by noting that when B is full row rank B T M~ l B is 
symmetric and positive definite. 

In our later controller design and stability analysis, we only use the fact that such 
constants Cmin and Cmu exist, while the actual estimation of these constants is not needed. 
To close this section, we re-state the control objective as follows: find a smooth control u 
without using the detailed information of P(.) and L(.) such that the system tracking error 
vector z locks in a compact set containing origin as time goes by, i.e., |]z|| < // as t > T 
(ji is a small constant related to control precision). 


3. Neuro-adaptive Fault-tolerant Control 

The nonlinear and uncertain nature of P(.) and L(.) represents the major challenge in 
control system design. The typical approach is to re-organize L(.) as a linear combination 
of some nonlinear functions. This has been the basis for most stable adaptive controllers, 
known as the regressor based method in literature [6]. For complex systems such as 
composite structures, this method involves complicated design procedures and heavy 
computations. In the early work of neural network based control, a common practice is 
to construct multilayer NNs and train these networks via gradient algorithms [7]-[9], 
[12]-[13] to approximate P(.) and L(.) respectively. Recognizing the potential instability 
associated with these algorithms, several researchers have proposed NN control schemes 
using on-line and stable training mechanisms based on Lyapunov stability theory [11]- 
[ 12 ]. 

In this work we develop a neuro-adaptive control scheme for the system (1). As a first 
step, we “reconstruct” L(.) via the first NN unit as 
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(7) 


Lnn = U) +Wi 7 Vi(z) 

where Lo represents the priori (crude) estimate of L, W\ e /? mxin is the “optimal” weights 
of the NN, t//i E /? m is the basis function whose selections typically include linear 
function, sigmoidal function, hyperbolic tangent function, radial basis functions etc. [7]- 

[15]. 

Let the discrepancy between L and L NN be denoted by 

£ = L-L NN (8) 

which is termed as NN reconstruction error in literature. With the support of the 
approximate theory of multilayer neural networks[7]-[8], such an error can be made 
arbitrarily small provided that certain conditions (e.g., sufficiently large number of 
neurons, sufficient smoothness of L etc.) are satisfied. This has been the basis for many 
NN based control laws [16]-[19]. 

Note that, however, in reality one can only construct a network with a finite number of 
neurons and that the nonlinear function being reconstructed is not smooth enough, as in 
the case of a jump fault. Thus the approximation capabilities of the network can be 
guaranteed only on a subset of the entire system state space. In other words, during the 
system operation, the reconstruction error is bounded by a constant sometimes, but 
unbounded at other times. For this reason, the NN reconstruction error should be treated 
carefully in control design because simply assuming that ||£j| is bounded by a constant 
does not warrant the success of the control scheme. For executing the idea, in our 
development, we remodel reconstruction error as follows: 

e = ce l +(l-c)e 2 (9) 


with 



if the NN works satisfactorily 
if the NN fails to work well 


To make the problem solvable, we make the assumption that 

II e \ IN < °°> 311(1 Ite INI ^-^-^iViIMiC). 
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Remark 2 


The above model seems to be more reasonable to reflect the fact that in general the 
reconstruction error could take any form, depending on many factors such as the basis 
function, the number of neurons, the structure of the networks, the weight updating algorithms 
as well as the original nonlinear function to be approximated. 

The point to be made is that do and d\{.) are still unavailable and should not be included in 
the control scheme directly. Another major challenge stems from the fact that we do not 
know when the NN behaves well. Sanner and Slotine [16] proposed a method of 
partitioning the state space into Ad and ^ and assumed that the NN works satisfactorily 
within Ad. The problem with this method is that it may be practically hard to find such 
subset Aj at which the NN works fine during the system operation. 

Here we take a different method to address this issue. Namely, we incorporate two NN 
units in the control scheme, with the first NN unit compensating the lumped 
nonlinearities/uncertainties and the second NN unit attenuating the effects due to the NN 
reconstruction error and the other resulting uncertainties. The overall control scheme is 
given by 

u = -kQZ + L NN +u c (10a) 

where 

Lr N =[^+w/m(z) (10b) 

denotes the control action of the first NN unit and u c is a compensating signal from the 
second NN unit to be specified. Here ko> 0 is a constant chosen by the designer, W is the 
estimation of the ideal weight values, Lq is the prior (crude) estimation of L and yf \ is 

again the basis function vector. With the control of the above scheme, the original 
system (6a) becomes 

z = -IcqPz+w/^i-Pu c +7J ( 11 ) 

~ A 

where Wj = Wj - Wj is weight estimation error and 
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( 12 ) 


rj = e+(I-P)(L 0 +W l T Vi) 

is the resultant uncertainty to be attenuated by the second NN unit. 

The design of the second NN unit consists of two steps. The first step is to formulate the 
“magnitude” of tj, which allows the m-dimensional nonlinear function to be “funneled” 
into a scalar nonlinear function. The second step is to build a NN to fight against the 
resulting scalar quantity, i.e., the “funneling effect” of 77 . For this reason the second NN 
is named as a funneling NN.We use the Euclidean norm to funnel the effect of tj, i.e. 

1^1 <y (13a) 

To find the expression for y, we make the assumption that there exists some unknown 
constants and some scalar nonlinear functions Y\ such that 

1=1 

where n, is an integer, which, together with (8) and (12), leads to 

N \+ 3 

y= E (oft (13b) 

i = 1 

with 

flty+i = d 0> <% 1+ 2 ®N l+ 3= SU P IU-^II. 

z sR c <zR m 

- *• y n,*2 =11 n II. s + 3 HI wm II 

Note that while Y, are known scalar functions, o\ are unknown in general, which makes y 
unobtainable. With the funneling NN, y is estimated as 

f=E^n(l/) = ^Vn(n ( 14 > 

i=l 

where N = N\ +3, d>, denotes the estimated weights and W\ui-) are some basis functions 
satisfying 


- 8 - 


Vm( Y i)- Y i 


(15) 


Based upon y , the control action of the second NN unit is given by 


u 


c 


X^Ui^VnOO 

i=l 

x ¥n x N +v 

i=l 


(16) 


where v > 0 is a small number chosen by the designer. To synthesize the tuning 
algorithms for Wj and (b , we consider the following Lyapunov candidate function 


v=v l +v 2 +v 3 


(17a) 


with 


v 'l=^ 7 '2. V 2 =^~,r(W l T W 0 , 

2 2 gl 

V 3 = (M-Cmin®) 7 (a - c min a)) 

^^II c min 


(17b) 


where gi > 0 and gn > 0 are free design parameters affecting weight learning rate and 
control performance, c,ni n > 0 is the minimum eigenvalue of P. Its derivative along the 
vector field ( 11 ) is 


with 


V = -Icqz T Pz + z^WjVi + z T r] - z T Pu c + V 2 + V 3 
^-VminNf+G +tr[w l T ( V , l z T -go'tii)] 


(18) 


Q = z t tj- z t Pu c — — (a ) - c min d>) T d) 

Sn 

Before presenting the weight tuning algorithms, let us focus on Q. With u c as in (16) and 
the upper norm bound on 7 as in (13a), we have 
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( 19 ) 


QAUWY 


c minS^ ni lUII ^Vn , 

i=l 1_ 

N a„ 

v i ir„ ii *r ii j.„ «n 


{(o-c^ n &) T & 


^We { lUII+v 

i=i 


= (co-c min d)) 


I VEi II Z II 2 Ve . 

i=l 1 

j 1 A1 


/V 


I^nj lUII+v 
L i=i 1 


(O 

£n 


+ V0), 


m 


where (Dm- max (ft*). Note that in deriving (19) we have used the facts that 


0) T Y < (J We <(0 m Y, We, 
i=l 


and 


2 Vu, II z || 

• t 1 


i=l 


N 


<1, Vv>0 


I^Hj IUII +V 

i=l 


Combining (18) and (19), we come up with the following weight tuning algorithms: 


W l =-a l W l + g x W\Z l 

N 


I II z || We 


d> = -a u d>+g u i=L 


( 20 ) 


I l^IIj II Z || +v 
i=l 


Consequently, 


^ - -c min^0 II 7 - |P +— (r(lV[ T lV|) 

Si 


+^(cd- c min a>) T & + vo), 
8n 


( 21 ) 


m 


Completing the square we can further rewrite (21) as: 
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( 22 ) 


V^-c mi „*olUH 2 -f l -<r(W I T lj' I ) 

2 S I 

- - — - - c min ") 7 - Crnin^) + * 

2 gn c min 

= -W-W2-ty3 + # 


where 


5 = vo) m +-^-tr^W l )+-^0) T a) 
2gi 2g n 


A = 2c min V ^ - °I, ^3 ~ °II 


(23) 


are constants independent of z. Note that (22) implies that ||z||,Wj and Q) are bounded. 
Furthermore, from (22), we have 


V<-fc 0 c min |zf+<y 

It is seen that v is negative as long as 


H> 


'_s ' 

V ^O c nun J 


* 


Therefore, z is confined in the region Q. = 


(24) 


z \z ^ 


.Vi 


% 


mm J 


(25) 


. From previous lemma. 


we establish that the control error is bounded. 


4. Simulation Verification 

Simulation on a numerical example is conducted to test the effectiveness of the proposed 
method. The parameters involved in the algorithms are chosen as 

*0=10, a i = (T/j =0.5, gj =gjj =0.5, v = 0.01, AT = 0.0 1 sec 
The basis functions for the first NN unit are chosen as: 
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yfi{c u ,X,a) 


l+e^Hi) 2 


(cn = 0.2, 0.4,..., (Xl = 


) 


where a = 0.5. The basis functions for the second NN unit are selected as in the 
following: 

\2 


yr n {c 2i ,X,a) = 


c 2i~ c 2i g 


-a(lX\-c2if 


/j + e 


-a{\x\-C2i) 2 


(c 2i = 0.3, 0.7 , ..., |XI = 


) 


where fi = 0.5. Ten NN basis functions were used for each of the NN 


units: iffeR J , j = 10 . All the weights are initialized at zero in the simulation. 

As mentioned earlier, any fault due to damage, element malfunction, and so forth, plays 
an important role in the behavior of the system. To have a realistic simulation it is 
necessary to model the aerodynamic effects of the damage. For this purpose, it may be 
said that the changed aircraft dynamics following damages to a stabilator surface is 
mainly due to an instantaneous change of the normal force coefficient of the surface, with 
the axial force coefficient being usually negligible and with the moment coefficient being 
proportional, through the geometry, to the normal force coefficient. Therefore we 
consider the expressions for fault in the simulation as follows. 


M f (.) = j3 0.3 M 0 
N f (.) = 0-( N { r + N 2 r + N 3 ) 

where /? = is the fault function, Th-15 sec is the time instant at which the 

fault occurs, m=0.2 is the fault coefficient. The results are presented in Figures 1-4, 
where Figure 1 is the stabilization error for critical and non-critical modes. Figure 2 
shows the control signals. Figure 3 and Figure 4 illustrated the on-line updating of 
weights for the two NN units. 
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Trajectory error for critical mode 



Figure 1 - System Performance Under Faults 


Control response 



Figure 2 - Control Input Under Faults 
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The plot of W 1 hat 



Figure 3 - Weights Updating for 1 st Neural Network 

The plot of W2hat 



Figure 4 - Weights Updating for 2 nd Neural network 


Also, for test purpose, we conducted a simulation on one critical mode with both jump 
and ramp faults. The results are presented in Figures 5-8, where Figure 5 is the control 
error under the ramp and jump faults occurring at 7y = 20 seconds. Figure 6 shows the 

control inputs and Figures 7 and Figure 8 are the updated weights of the networks. As 


! 

i 
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can be seen, with the proposed control method the system remains stable under severe 
faults and all the internal signals are bounded and smooth. 




a) ramp fault b) jump fault 

Figure 5 - System Performance Under Faults 





a a to IB 2< 



Thn«(*) 


b) jump fault 
Input Under Faults 


CoaOcMrt* of F Mt Nmuml H mtwotk 



7 o ao «o 100 


fer Ihi Second NmmI 



a) ramp fault b) jump fault 

Figure 7 - Weights Updating for the 1 st Network I 
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a) ramp fault b) jump fault 

Figure 8-Weights Updating for the 2 nd Network II 


5. Conclusion and Future Work 

Whenever an actuator fault occurs, the element of B is of the form, 

9t ..= ! 1 

where 




0 

- ( 0 , 1 ) 

1 


|<r ,70 ± fi(t -Tf )a u (.)| = 0 implies that the actuator is totally failing 
j <T t ,0 i P(t-Tf )<t, 7 (.)| = 1 implies the actuation system is health 

|<r ,70 ±($(t-Tf )or,7 (.)| = S < 1 the fading actuation system is involved 
We are currently investigating the control problem under the following fault conditions, 

0 < |otfo ±P(t-T f )a u (.)| < 1 

This, physically, involves a fading actuation fault in which the actuating system has not 
totally failed yet. The results will be reported in the forth coming publication. 
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